Robust inverse energy cascade and turbulence structure in three-dimensional layers 
of fluid 

D. Byrne, 1 H. Xia^Qand M. Shats 1 

Research School of Physics and Engineering, 

The Australian National University, Canberra ACT 0200, 

Australia 

(Dated: 14 September 2011) 

Here we report the first evidence of the inverse energy cascade in a flow dominated 
by 3D motions. Experiments are performed in thick fluid layers where turbulence 
is driven electromagnetically. It is shown that if the free surface of the layer is not 
perturbed, the top part of the layer behaves as quasi-2D and supports the inverse 
energy cascade regardless of the layer thickness. Well below the surface the cascade 
survives even in the presence of strong 3D eddies developing when the layer depth 
exceeds half the forcing scale. In a bounded flow at low bottom dissipation, the 
inverse energy cascade leads to the generation of a spectral condensate below the 
free surface. Such coherent flow can destroy 3D eddies in the bulk of the layer and 
enforce the flow planarity over the entire layer thickness. 



I. INTRODUCTION 

There has been remarkable progress in the 
understanding of turbulence in fluid layers. 
Such layers, characterized by large aspect ra- 
tios, are ubiquitous in nature. 2D turbulence 
theory by Kraichnan 1 , in particular the in- 
verse energy cascade, has been confirmed in 
experiments in thin fluid layers^l More re- 
cent studies showed that theoretical results 
derived for idealized 2D turbulence are valid 
in a variety of conditions, even when the the- 
ory assumptions are violated^. 

In bounded turbulence, the inverse en- 
ergy cascade may lead to the accumulation of 
spectral energy at the box size scale and the 
generation of a spectral condensate, a large 
vortex coherent over the flow domairP. Good 
agreement with the Kolmogorov-Kraichnan 
theory was found in the double lay ers o f fluids 
in spectrally condensed turbulenc e 1 9 1 10 1 . 

It should be noted that though spectral 
condensation was ob serv ed in the double- 
layer configur at ions^^ 9 ^, in single layers of 
electrolytes not only spectral condensation, 
but also the very existence of the inverse en- 
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ergy cascade has been questioned. For ex- 
ample, in Ref.^ the flow generated in a sin- 
gle layer was referred to as "spatio-temporal 
chaos" to stress the absence of the turbu- 
lent energy cascades. It is thus important to 
better understand spatial structure of turbu- 
lence in such layers as well as differences be- 
tween turbulence in double and single layers. 
(Here we do not discuss MHD flows which 
also show spectral condensation^, but where 
two-dimensionality is imposed by homoge- 
neous magnetic field and bottom dissipation 
is restricted to a thin Hartmann layer.) 

A comparison of turbulent flows in a single 
and double-layer configuration is also impor- 
tant to improve our understanding of turbu- 
lence in atmospheric boundary layers. These 
layers are very different over terrain and the 
oceans, with the former being substantially 
thicker than the latter ones. Stable immis- 
cible layers of fluids have been generated 
in the laboratory by placing a heavier non- 
conducting fluid at the bottom of the cell and 
a lighter layer of electrolyte resting on top of 
fiPHnD In this case the electromagnetic forc- 
ing is detached from the solid bottom and it 
is maximal just above the interface between 
the two fluids. The structure of the flow in 
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the top layer close to the interface with the 
bottom layer may be similar to that of the 
atmospheric boundary layer over the ocean. 

In this paper we present new results on 
the spatial structure of turbulence in a single- 
and a double-layer turbulent flow. In con- 
trast to the previous studies, we focus here 
on thick layers. Recent 3D direct numerical 
simulations of the Navier-Stokes equations^ 
have shown that the finite layer depth leads 
to splitting of the energy flux. In thin lay- 
ers the energy flux injected by forcing is in- 
verse. As the layer thickness h is increased, a 
larger fraction of the flux is redirected down 
scale, towards the wave numbers larger than 
the forcing wave number kf. At h/lf > 0.5, 
most of the injected flux cascades to small 
scales. 

Experiments in thick single layers of elec- 
trolyte partially confirm this picture 1 - 4 . It has 
been shown that at h/lf > 0.5 strong 3D mo- 
tions indeed appear in a layer. Such 3D ed- 
dies increase vertical flux of the horizontal 
momentum (eddy viscosity) which leads to a 
bottom drag dissipation rate higher than the 
one expected from a quasi-2D model. How- 
ever, no evidence of the direct energy cascade 
was found. At k > kf no fc~ 5//3 spectrum was 
observed as should be expected in the pres- 
ence of the direct energy cascade. Probably 
3D eddies over the solid bottom introduce 
higher damping for small scales (k > kf) and 
do not allow the direct cascade to develop. 

In the double layers the effects of the solid 
bottom on the top layer turbulence are iso- 
lated by a layer of non-conducting, low vis- 
cosity heavier fluid at the bottom. It was 
recently founcP^ that the range of depths in 
which a top layer flow remains planar and 
supports spectral condensation, is greatly ex- 
tended, h/lf > 0.5. It was hypothesized, that 
a residual inverse energy flux condenses tur- 
bulent energy into large scale flows even in 
thick layers. These coherent flows enforce the 
planarity by shearing off vertical eddies and 
thus secure the upscale energy transfer^. It 
is not clear however how the large scale flow 



is generated in a thick layer in the first place. 
In this paper we present new results which 
shed light on this. 

II. EXPERIMENTAL SETUP AND 
FLOW CHARACTERIZATION 

In these experiments turbulence is gener- 
ated through the interaction of a large num- 
ber of electromagnetically driven vortices^'. 
An electric current flowing across the con- 
ducting fluid layer interacts with the spatially 
varying vertical magnetic field produced by 
a 24 x 24 or 30 x 30 array of magnetic 
dipoles (10 mm and 8 mm separation respec- 
tively), Figure [TJ The magnet arrays are 
placed under the bottom of the 0.3 x 0.3 m 2 
fluid cell. To ensure that turbulence is forced 
monochromatically at k = kf, vertical mag- 
netic field produced by the array has been 
measured using a Hall probe scanned in hori- 
zontal planes at several heights above the ar- 
ray. The measured magnetic field (Fig. [ljb)) 
has then been Fourier transformed in 2D, 
Fig. [l](c) . The spectrum shows that J x B 
forcing is localized in fc-space in a narrow 
spectral range (in this example, for a 10 mm 
magnet separation, the spectrum peaks at 
k pa 630 rad-m -1 ). The forcing strength 
is controlled by adjusting electric current 
through the layer in the range 0.5-5 A. 

Turbulence is generated either in a single 
layer of Na2SO^ water solution, or it is driven 
in the top layer of the electrolyte which rests 
upon a layer of heavier (1820 kg/m 3 ) non- 
conducting liquid (FC-3283 by 3M) which is 
not soluble in water. In the latter case, forc- 
ing is the strongest just above the interface 
between the layers. 

The flow is visualized using neutrally 
buoyant seeding particles 50 /zm in diameter 
illuminated by a horizontal laser sheet. The 
thickness of the laser sheet and its height rel- 
ative to the free surface are adjusted to vi- 
sualize different regions of the layers. Per- 
turbations to the fluid surface are monitored 
by reflecting a thin laser beam off the free 
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FIG. 1. Schematic of experimental setup, (a) Neutrally buoyant seeding particles in the top 
(conducting) layer are illuminated and their motion is filmed from above, (b) Measured vertical 
magnetic field produced by the magnet array: blue and red dots indicate upward and downward 
B direction, (c) Wave number spectrum of the measured magnetic field. 



surface onto a distant screen (the detection 
sensitivity is ~ 5 x 1CT 3 mm). In all reported 
experiments no perturbation to the free sur- 
face is detected. 

To characterize the vertical structure of 
the flow, defocusing PIV is usecP^. This 
technique uses a single CCD camera with 
a multiple pinhole mask to measure three- 
dimensional velocity components of individ- 
ual seeding particles in the flow^. The hor- 
izontal turbulent velocity fields are derived 
using particle image velocimetry (PIV). Par- 
ticle pairs are matched from frame to frame 
throughout the illuminated volume using a 
PIV/PTV hybrid algorithm. Derived ve- 
locities are then averaged over hundreds of 
frame pairs to generate converged statistics 
of the root-mean-square velocities (V Xyy>z ) 
throughout the layer. The technique allows 
one to resolve vertical velocities (V z ) rms > 0.4 
mm/s and (V xy ) rms > 2.5 xlO -2 mm/s. The 
damping rate of the flow is studied from the 



decay of the horizontal energy density after 
switching off the forcing^!. 

We also study vertical motion of seeding 
particles by illuminating the layer using a 
thick vertical laser sheet. Streaks of the par- 
ticles in z — x plane are filmed with the expo- 
sure time of 1-2 s through a transparent side 
wall of the fluid cell. 

III. TURBULENT FLOW IN A 
THICK SINGLE LAYER 

We first describe turbulent flows in a sin- 
gle layer. The flow is forced at kf rs 800 
rad-m -1 (// ~ 7.8 mm) in a layer of thick- 
ness hi = 10 mm. According to numerical 
simulations, turbulence in such a layer should 
show substantial three-dimensionality^, even 
when the forcing is 2D. Figure [2] shows verti- 
cal profiles of vertical (a) and horizontal (b) 
velocities along with the snapshot of the par- 
ticle streaks in the vertical (z — x) plane. 
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FIG. 2. Vertical profiles of (a) vertical, V z , and 
(b) horizontal, V XiV , velocities. A grey box in (a) 
indicates the sensitivity of the defocusing PIV 
technique, (c) A snapshot of the particle streaks 
taken at the exposure time of 2 s. 

RMS vertical velocities < V z > are low in the 
top sublayer, 2 mm below the free surface, as 
well as in the bottom boundary sublayer. In 
the bulk of the flow (2-8 mm) < V z > is 
high, being only a factor of two lower than 
horizontal velocities < V x>y >. 

< Vx,y > shows a maximum at h — (3 
- 4) mm, which is indicative of the com- 
petition between the forcing and the bot- 
tom drag. The forcing / ~ (J x B) is the 
strongest near the bottom (magnets under- 
neath the cell, the current density J is con- 
stant through the layer) and it decreases in- 
versely proportional to the distance from the 
bottom, / oc 1/h. The bottom drag is also 
the strongest near the bottom; in a quasi-2D 
flow it should scale faster^, a oc 1/h 2 , result- 
ing in the maximum of < V x>y > at h — (3—4) 
mm. In the top sublayer, h% = (8 — 10) mm, 
turbulence is expected to behave as quasi-2D 
due to the lower vertical velocities and the 
absence of vertical gradients of the horizon- 
tal velocities. The planarity of the top sub- 
layer can be seen qualitatively in the particle 
streaks of Fig. W[c). 



To test if the nature of the turbulent en- 
ergy transfer changes between the top sub- 
layer, the bulk flow and the bottom layer, we 
perform PIV measurements of the horizon- 
tal velocities by illuminating three different 
ranges of heights: hi — (8 — 10) mm, h2 = 
(4 — 7) mm, and h2 — (1 — 4) mm. The ve- 
locity fields are analyzed as described, for ex- 
ample, in RefP. We compute the wave num- 
ber energy spectra Ek(k) and the third-order 
structure function S 3 l = ((SVl) 3 ), where SVl 
is the increment across the distance r of the 
velocity component parallel to r. The third- 
order structure function is related to the en- 
ergy flux in fc-space as e = — (2/3) S^l/?"- Pos- 
itive S%l corresponds to negative e and the 
inverse energy cascade. 

Figures (J3|(a,b) show the kinetic energy 
spectrum and the third-order structure func- 
tion function of the separation distance 
measured in the top sublayer h\. At k < kf 
the spectrum scales close to k~ 5 ^ 3 , while S^l 
is positive at I > If and is a linear function 
of r. This is in agreement with the expec- 
tation of the quasi-2D turbulence in the top 
sublayer. 

In the bulk flow h2, which is dominated 
by 3D motions, the spectrum is still close to 
k~ 5 ^ 3 , though it flattens at low wave num- 
bers as seen in Figj3^c). Consistently with 
this, the range of scales for which S$l is pos- 
itive and linear is reduced to about r w 40 
mm, Fig|3|d). Such a behavior of S3z(r) is 
also typical for thin single layers. The re- 
duction in the inverse energy cascade range 
is correlated with the increased damping. 3D 
motions present in the bulk flow (Fig|3^c)) in- 
crease damping rate due to the increased flux 
of the horizontal momentum to the bottom of 
the eel™. The increased damping arrests the 
inverse energy cascade at some scale smaller 
that the box size. 

In the bottom sublayer h$, the flow is sub- 
ject to even stronger damping. As a result, 
the spectrum is much flatter than k~ 5 ^ 3 . S3L 
is positive for a narrow range of scales, giving 
a hint of same trend as in the bulk and the 
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FIG. 3. (a,c,e) Wave-number spectra and (b,d,f) the third order structure functions S^l measured 
in (a,b) the surface sublayer (h\ = (8 — 10) mm), (c,d) in the bulk flow (/12 = (4 — 7) mm), and 
(e,f) in the bottom sublayer regions (/13 = (1 — 4) mm). 



top sublayers. 

The above results suggest that despite the 
presence of substantial 3D motion in a thick 
(h/lf = 1.28) single layer, statistics of the 
horizontal velocity fluctuations remain con- 
sistent with that of quasi-2D turbulence and 
supports the inverse energy cascade. As seen 
from the energy spectra of Fig|3^a,c,e), there 
is no evidence of the direct energy cascade at 
k > kf. The spectrum shows that scales 
much steeper than the 3D Kolmogorov scal- 
ing of &;~ 5//3 . A possible reason for the ab- 
sence of the direct energy cascade in a 3D flow 
at k > kf is the fact that the Reynolds num- 
ber is low in these experiments (Re < 300). 

We performed experiments in even thicker 



layers, up to h/lf rs 2.3. If the layer thickness 
is further increased in a bounded flow, a spec- 
tral condensate forms in the top sublayer h\. 
In a layer of a total thickness of 20 mm the 
spectral condensate penetrates down to 4-5 
mm below the free surface. The formation of 
the condensate in thick layers is due to the 
reduction in the bottom damping, since even 
in the presence of 3D motions the damping 
rate is reduced with the increase in thickness 
as 231 a ~ l/h. 
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FIG. 4. Vertical profiles of (a) vertical, and (b) horizontal RMS velocities measured at t = 5 s 
after forcing start, solid diamonds, and at t = 20 s, open circles. 



IV. TURBULENCE STRUCTURE 
IN A DOUBLE LAYER FLOW 

As shown in the previous section the in- 
verse energy cascade is sustained in the pres- 
ence of 3D motions. In bounded turbulence 
at low damping such a cascade leads to spec- 
tral condensation and to the generation of a 
large-scale coherent structure. Such a struc- 
ture could then impose two-dimensionality on 
the flow in the layer, as has been found in re- 
cent experiments^. In this section we study 
spatial structure of the flow during spectral 
condensation in thick layers subject to even 
lower bottom drag. 

The bottom drag is reduced by generat- 
ing two immiscible layers in which the bot- 
tom layer (heavier, non-conducting liquid) 
isolates the conduction layer from the bot- 
tom. We keep the bottom layer relatively 
thin, hb/lf < 0.5. The top layer, on the other 
hand, is thick, ht/lf > 0.5, to allow three- 
dimensionality to develop, even with 2D forc- 
ing. Here we study the layer configuration 
described in Ref.™ namely h t = 7 mm and 
hi — 4 mm, which correspond to hb/lf = 0.44 
and ht/lf = 0.78 respectively. It has been re- 
ported that the flow in the top layer shows 
substantial 3D motions shortly after turbu- 
lence is forced. However in the steady state, 
the development of spectral condensate leads 
to a substantial reduction of 3D eddies and 
to a planarization of the flow. 



Figure [4] shows vertical profiles of (a) 
vertical, and (b) horizontal RMS velocities 
measured using defocusing PIV in the top 
layer. The grey box in FigQa) indicates the 
method sensitivity. Shortly after the flow is 
forced (t =5 s), vertical velocities peak in the 
middle of the layer at (V z ) ~ 1.7 mm/s, solid 
diamonds in Fig|l](a). 20 seconds later verti- 
cal velocity fluctuations are substantially re- 
duced, to below the resolution level, (V z ) < 
0.5 mm/s. Horizontal velocities during the 
initial stage of the flow development peak at 
the interface between the two layers, i.e. in 
the region of strongest forcing. In the steady 
state, after the large coherent vortex develops 
the horizontal velocity profile becomes flat 
over most of the top layer, showing (V X)V ) ~ 
6 mm/s at h = 2 — 7 mm above the inter- 
face, open circles in FigQb). A reduction in 
(y x , y ) near the interface in the steady state 
is probably related to the effect of sweeping 
of the forcing scale vortices by the developing 
condensate^. Thus, a substantial fraction of 
the top layer is quasi-2D, i.e. (V z ) ~ and 
d(V XtV )/dz*i0. 

Such a quasi-two-dimensionality of the 
flow is attributed to the shearing of verti- 
cal eddies by a strong condensate^. Indeed, 
in the steady state, the statistics of horizon- 
tal velocities is in a good agreement with the 
Kraichnan theory. The spectra and the third- 
order structure functions in the presence of 
spectral condensate are computed after sub- 
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FIG. 5. Statistics of horizontal velocities in a 
double layer configuration: top layer thickness 
ht = 7 mm, bottom layer thickness ht = 4 mm. 
The forcing scale is If = 7.8 mm. (a) Wave 
number spectrum of the kinetic velocity Ek, and 
(b) the third-order structure function S^l, com- 
puted after subtracting time-averaged mean ve- 
locity field. The entire top layer is illuminated. 

tracting the time-average velocity field from 
the instantaneous velocity field, as discussed 



jrPtLp] After the mean subtraction, the Spee- 
ds 



trum shows Ek ~ /c _5/s , as seen in Figj5[a) 
The third-order structure function is positive 
and is a linear function of the separation dis- 
tance r up to r 70 mm, Figjf^b). Thus the 
spectra and the structure functions in such 
a flow agree with quasi-2D expectations and 
are consistent with the vertical structure of 
the flow of Fig. |4| 

V. SUMMARY 

We have studied spatial structure of tur- 
bulent flows in thick layers at low Reynolds 
numbers {Re ~ 100). If the free surface of a 
layer is unperturbed, there is a finite thick- 
ness layer close to the surface, which remains 
quasi-2D regardless of the total layer thick- 
ness. Two-dimensional turbulent velocities 
in the top layer show kinetic energy spectra 
and the third-order structure functions con- 
sistent with the Kraichnan theory of 2D tur- 
bulence. Vertical eddies which appear when 
the layer is sufficiently thick, h/lf > 0.5, in- 
troduce additional bottom drag due to the 
eddy viscosity^, but they do not qualita- 



FIG. 6. Schematics of the structure of turbu- 
lence in thick layers. 

tively change the statistics of the horizontal 
velocity fluctuations, which remains quasi-2D 
even in the presence of 3D motions. If the 
bottom drag is reduced by introducing an im- 
miscible thin bottom layer, the inverse energy 
cascade leads to spectral condensation and 
to the formation of the large scale coherent 
structures. Such flows, as has recently been 
showr/M shear off eddies in the vertical plane 
and reinforce quasi-two-dimensionality of the 
flow. Measurements presented here, in par- 
ticular Fig|4j confirm this. 

Summarizing, flows in thick layers of flu- 
ids with an unperturbed free surface can be 
viewed as two interacting sublayers, as illus- 
trated in FigJHj The top layer is quasi-2D; 
it supports the inverse energy cascade. In 
a bounded domain at low damping, the in- 
verse cascade leads to spectral condensation 
of turbulence. The bottom sublayer is dom- 
inated by 3D motions which are responsible 
for the onset of the eddy viscosity. A planar 
coherent flow (spectral condensate) develop- 
ing in the top layer can reduce the bottom 
layer thickness /12 through shearing of the 3D 
eddies. The thickness of the two sublayers 
thus depends on the competition between the 
vertical shear and the 3D motions due to the 
forcing. In the two layer configuration, the 
spectral condensate formed in the top sub- 
layer can take over almost the entire layer 
thickness. 
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